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PREFACE 


This report outlines a procedure for estimating the ultimate adjustment of 
the shoreline and of the offshore bathymetry to changes in mean water level on 
the Great Lakes. The procedure is based on a sediment balance model, cali- 
brated and verified using profile changes measured on the eastern shore of 
Lake Michigan over a 9-year period. The procedure is generalized for applica- 
tion to other sections of the Great Lakes by considering regional variations 
in storminess, and by requiring local evaluation of relevant geomorphic and 
textural variables. 


The report was prepared by Edward B. Hands, under the supervision of 
Dr. C.He Everts, Chief, Engineering Geology Branch. Dr. W.Le Wood of Purdue 
University and J. Pope of the U.S. Army Engineer District, Buffalo, provided 
data used in the example problems. Appreciation is also extended to 
We Birkemeier, C.H. Everts, R.J. Hallermeier, R.D. Hobson, Re Jachowski, and 
E. Meisburger for their review and suggestions which significantly improved 
this report. 


Comments on this publication are invited. 


Approved for publication in accordance with Public Law 166, 79th Congress, 
approved 31 July 1945, as supplementad by Public Law 172, 88th Congress, 
approved 7 November 1963. 


ED Ee BISHOP 


Colonel, Corps of Engineers 
Commander and Director 
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CONVERSION FACTORS, U.S. CUSTOMARY TO METRIC (SI) UNITS OF MEASUREMENT 


U.S. customary units of measurement wee in this report can be conwentcll to 
metric (SI) units as follows: 


Multiply 


To obtain 


ES 


inches 


square inches 
cubic inches 


feet 


square feet 
cubic feet 


yards 
square yards 


cubic yards 


miles 
square miles 


knots 

acres 
foot-pounds 
millibars 
ounces 


pounds 


ton, long 
ton, short 


degrees (angle) 


Fahrenheit degrees 


25.4 
2.54 


453.6 


& 10°82 


millimeters 
centimeters 
square centimeters 
cubic centimeters 
centimeters 
meters 
square meters 
cubic meters 
meters 
square meters 
cubic meters 


kilometers 
hectares 


kilometers per hour 

hectares 

newton meters 

kilograms per square centimeter 
grams 


grams 
kilograms 


metric tons 
metric tons 
radians 


Celsius degrees or Kelvins! 


Te 


lTo obtain Celsius (C) temperature readings from Fahrenheit (F) readings, 


use formula: 


To obtain Kelvin (K) readings, 


(5/9) (F -32). 


use formula: 


(S/9) (CF Ss) se 2WSo NS 


SYMBOLS AND DEFINITIONS 
sediment overfill factor--the ratio of sediment volume supplied by 
profile recession to that retained after sediment sorting, packing, 


and profile readjustment 


signum function having values of: 1 for z > 0; =l for z < 0; and 0 
for z = 0 


average horizontal extent of profile adjustment 
average horizontal displacement of the profile 
average vertical extent of adjusting shore profile 


average change in elevation of the bottom profile or the water surface 


PREDICTING ADJUSTMENTS IN SHORE 
AND OFFSHORE SAND PROFILES ON THE GREAT LAKES 


by 
Edward B. Hands 


I. INTRODUCTION 


This report briefly describes a method for predicting long-term changes in 
shoreline position and offshore bathymetry on the Great Lakes. Beach profiles 
fluctuate in response to storms and water level changes. On the Great Lakes 
both storms and water levels undergo prominent seasonal fluctuations. Super- 
imposed on the seasonal fluctuations is a longer term variation in annual mean 
water level elevations. When measured over a number of years the net long- 
term change in water levels exceeds the range of seasonal fluctuations. The 
method described here for predicting long-term profile adjustments to changing 
lake levels is based on a conceptually sound, empirically verified model which 
includes allowances for regional variations in storm exposure, coastal geomor- 
phology, and sediment texture. 


Il. THE IDEALIZED MODEL 


As described by Bruun (1962), a rise in the mean elevation of the water 
surface tends to shift the equilibrium sand profile landward. As water levels 
rise, erosion prevails on the upper beach, and the shoreline retreats. Con- 
ceptually, the erosion supplies material to build the outer part of the 
responding profile upward. Eventually, the initial profile shape is reestab- 
lished farther inland and, at a distance above its initial position, equal to 
the change in water level, 2z, as depicted in Figure l. If there are no 
longshore losses the ultimate retreat of the profile x can be calculated 
given the dimensions of the responding profiles, X and Z, and a measure of 
the stability of the shore-eroded material, Ry. 


2X (Ra) °8 Z 
Z 


where sg z= 1 if z > O, and sg z = -l1 if z < 0. Hands (1980) provides a more 
detailed description of the sequence of profile changes leading back to equi- 
librium and the derivation of the equation as a direct consequence of the 
conservation of sediment volumes. Also provided is an extended equation to 
cover the case of longshore imbalances in sediment transport. Though the 
concept behind the equation is straightforward, its evaluation in the field is 
problematic because the required dimensions of the responding profile (X and 
Z) will usually be unknown and depend on the local wave climate. Monitoring 
of beach and offshore changes in Lake Michigan has both verified the per- 
tinence of the equation and simplified its evaluation for sandy shores 
throughout the Great Lakes. | 


III. EVALUATION OF TERMS IN THE EQUATION 
1. Change in Water Level, z. 


This is the given or independent variable; it refers to the change in mean 
elevation of the water surface which disturbed the equilibrium of the beach. 


Upper Point of Profile Adjustment 


ee Profile 
UP 


rofile Adjusted To Submergence, z 


Elevated Water Surface 


Initial Water Surface 


Point of Profile 
> Closure 


Figure 1. Sketch of profile measurements required to predict shore adjustment 
to a change in water level elevation. Providing there is no net 
gain’ or loss outside the control volume, constancy of profile shape 
requires that the ultimate shore retreat x be equal to 2X/Z. 


The change in water level elevation should persist at least long enough for 
equilibrium to be reestablished. The time required to do so will depend on 
the magnitude of the water level change and on the occurrence of storms which 
provide the energy necessary for reshaping the shore. To be consistent with 
the time frame in which other terms of the equation will be evaluated, the 
mean water level surface is expected to remain at its new elevation for more 
than 1 year. The new water surface elevation may refer to expected changes 
resulting from proposed modification to the lake control plans, or it may 
refer to a past change in lake level which could have resulted from natural 
variations in the water supply, but whose effect remains to be determined. 


2. Height of the Responding Profile, Z. 


This term refers to the vertical relief of the active beach which extends 
from the closure depth offshore to the top of the affected deposits on the 
backshoree The profile closure depth, beyond which the bottom does not 
respond to surface changes, depends on the local wave climate. Estimates of 
closure depths, based on the relationship between published wave climates and 
repetitive profiling on Lake Michigan (Hands, 1980), are tabulated in the 
Appendix for all five Great Lakes. The average height of the affected back- 
shore deposits above the initial water surface should be determined from field 
measurements at the particular site of application. The sum of the closure 
depth (taken from the App.) and the average backshore height (determined from 
field surveys) is an estimate of Z. 


3. Width of the Responding Profile, KX. 


The width of the responding profile refers to the horizontal distance 
between the landwardmost point of profile adjustment and the offshore closure 
depth. The width along many shore-normal profiles in the region of interest 
should be measured and averaged to obtain a representative value (Fig. 1). 


4. Overfill Ratio, Ry. 


The percent of eroded material by volume to be carried as suspended load 
beyond the closure depth should be estimated; additional erosion must compen- 
sate for the loss from the active profile. Hobson (1977) explains how to 
compute R,, based on textural parameters of “native beach” and “borrow” 
materials. The same procedures can be applied here except that the parameters 
for the borrow materials must be based on a composite sample of the eroding 
section of shore; ieee, the upper beach, in the case of an increase in lake 
level, because this is the zone which supplies sediment to rebuild the profile 
(Fige 1). If lake levels decline, then erosion below the lever point supplies 
material to prograde the upper part of the profile (Fig. 1). Im this case, 
the lower profile (between lever and closure points) corresponds to the 
“borrow area.” In either case “native” characteristics must be based on a 
composite sample of the entire responding profile, from the upper limit of 
profile adjustment to the point of profile closure. The lever point is useful 
in describing the sediment balance concept. In practice it is an ill-defined 
transition zone separating areas of predominant~ erosion and deposition. For 
long-term adjustments this transition may occur in the vicinity of the outer 
bar, but its exact location is not critical for the application of present 
procedures. 


If adequate data for calculating R, are not available and the lake level 
is rising, Ry, can be assigned a value of 1, provided that (a) there is less 
silt and clay in the beach and backshore than there is offshore, and (b) the 
mean grain size across the beach and backshore is greater than the mean size 
offshore. A value of Ry = 1 indicates that all eroded material is expected to 
remain in the zone of profile adjustment; if only P percent remains, then 
Ries 100/P. 


IV. EXAMPLE PROBLEMS 


The following problems are evaluated on the basis of limited available 
survey data. They provide examples of the basic steps in applying the pro- 
posed method of profile prediction. If these predictions were intended to 
support actual design or management decisions, a more careful evaluation of 
field conditions would be required. 


kk kk kk kk OK OK OK kK KOK &K EXAMPLE PROBLEM 1 * * * * & & & & & KOK OK KOR OK 


GIVEN: A contemplated change in the regulation plan controlling the water 
supply to Lake Ontario would raise the long-term surface elevation 0.3 
meter. 


FIND: The effect the higher stages would have at the eastern end of Lake 
Ontario. 


ANALYSIS: The barrier beaches and high dunes which characterize this stretch 
of shore are of special ecological and scenic value. Situated downwind from 
the major storm paths across Lake Ontario these barrier beaches are exposed 
to the highest storm waves reported on the Great Lakes, but because of 
relatively low land development there are few protective structures along 
this reach of the shore. Sand extends lakeward across a series of longshore 
bars. There are no known rock outcrops, and there is a close balance 
between southward and northward longshore transport. 


EVALUATION OF TERMS: 


z = 0.3 meter Proposed long-term increase in lake level. 


7-6 meters Average height of the eroding dunes above 
mean lake level (from field surveys). 


13.4 meters Profile closure depth (from the App-). 
Z = 21.0 meters Sum of the two values obtained above. 
X = 2,414 meters Average distance of the 13.4-meter depth 


contour from shoree The vertical datum 
should be the same as the reference level 
below which closure depth was measured in 
the previous step. 


Ry = il All of the material eroded from the upper 
beach is expected to remain within the 
bounds of the responding profile. 


It is estimated that the higher stages would shift the equilibrium shore 
profile an average of 34 meters inland and raise it 0.3 meter above its 
present elevation. 


kk kk k kK Kk KOK kK KK & & EXAMPLE PROBLEM 2 * * * & & & & kK KK KK KK 


GIVEN: Assume a new regulation plan is proposed to modify the inflow to Lake 
Michigan and Lake Huron via the St. Marys River. If adopted, this plan 
would lower the long-term mean surface elevation of Lake Michigan and Lake 
Huron by 0.3 meter. 


FIND: The effect the lower water levels will have on shore erosion at Indiana 
Dunes National Seashore. 


ANALYSIS: The dredged channel and navigation structures at Michigan City, 
updrift of the Indiana Dunes National Seashore, block some of the potential 
sediment input from the east. Westward longshore transport out of the dune 
area thus creates a sand deficit and contributes to a long-standing erosion 
problem in the park. As lake levels fall, the shoreline withdraws and the 
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beach widens. Assuming lake currents and waves are not altered, they will 
tend to reestablish the previous profile shape at a lower and more lakeward 
position. Longshore losses to the west continue to exceed the net supply from 
the east. However, offshore where the bottom slope is gradual, lowering of 
the water surface brings bottom sediments into a shallower hydraulic regime. 
This results in landward sediment transport which steepens the nearshore 
slope, builds dunes on the widened beach, and feeds the longshore currents 
leaving the dune area to the west. The cumulative effect of these adjustments 
can be estimated using the equation. 


EVALUATION OF TERMS: 


z = -0.3 meter Given (negative indicates a reduction in water 
level). 
2.9 meters Estimated average height of dunes expected to 


form on the widened beach lakeward of the 
present foredune. 


11.1 meters Profile closure (mean of depths at adjacent 
sites 28 and 29 in App-)- 


Z = 14.0 meters Sum of the two values obtained above. 


X = 3,030 meters Average distance of the 11.l-meter contour 
from shore, based on field surveys. 


Ra = i Offshore sands are expected to move, onshore, 
and the wind is not expected to carry sand 
inland past the present foredune. 


sg Zz Bs 
2X (Ry) 993,030) 1-4 
en ee = Oo meters (exalluaitimemthnemequattslon)) 


It is thus estimated that lowering the lake level 0.3 meter will effectively 
shift the equilibrium position 65 meters lakeward. As previously mentioned, 
a net loss of sand will still prevail due to the predominance of transport 
to the west. Therefore, the actual shoreline is not expected to advance 65 
meters lakeward. A reasonable interpretation is that there will be a long- 
term gain of 65 meters of beach that otherwise would have been lost by 
erosion if the water level had not been lowered. Dividing 65 meters by the 
average past recession rate would provide an estimate of when the avoided 
erosion would otherwise have, occurred. 


If 15 percent of the offshore sediments were thought to be too fine to 
remain in the active shore zone, then the width of shore “saved” should be 
reduced to (1 - 0.15) x -65 meters = -55 meters. Note that a liberal 
estimate of future dune heights would also make the predicted savings in 
beach width more conservative. 


kk kK kK Kk KKK KK KK KKK KKK KK KKK KK KK KK KK KK KK KK KK 
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Ve. LIMITATIONS 


The prediction scheme proposed here is based on a physically sound 
principle; its application has been shown to produce results in remarkably 
good agreement with actual measurements along a 50-kilometer reach of Lake 
Michigan's eastern shoree The study area included a wide range of shore types 
though sand was always the predominant material, and the range of wave con- 
ditions within the study area was narrow relative to its potential range 
among other sites of application (Hands, 1980). The model's simplicity 
should promote its widespread application. However, this very simplicity when 
contrasted to the actual complexities of nearshore processes, should also 
underscore the need for careful consideration of each application. Two broad 
areas of concern are (a) the interpretation of the results, and (b) the 
evaluation of the input values used in the prediction equation. 


1. Interpretation of Results. 


Recall that the model attempts to evaluate only the profile change induced 
by the change in water level. Other factors resulting in shore modification 
may be found in the Shore Protection Manual (U.S. Army, Corps of Engineers, 
Coastal Engineering Research Center, 19/77). If the other factors are signifi- 
cant, the resulting net effect can be combined with the water level induced 
change (as in eqe 4 in Hands, 1980). 


The key element in the proposed prediction scheme is the assumption that 
the profile will return to a specific shape after being disturbed by a change 
in water level. It is fundamental, therefore, that the reasonableness of this 
assumption be fully considered in each application. The proposed model does 
not apply, for example, to perched beaches, rocky shores, or beaches overlying 
nonerodible substrates that will be exposed during profile adjustment. 


To clarify the development of the equation, the initial and final profile 
shape has been referred to as the equilibrium profile. However, the known 
profile used to evaluate Z and X need not be in strict equilibrium so long 
as it represents a relatively stable shape to which the profile evolves after 
perturbation. It is not necessary for the profile to be stable before pertur- 
bation; e.g-, the shore in example problem 2 was continually eroding, both 
before and after the lake level disturbance. 


Obviously, seasonal and storm features also affect the profile shape. 
This complicates testing the assumption of constant overall shape, but it 
should not significantly affect the results of applying the discussed 
model. Of course, if significant seasonal fluctuations are known, they can be 
included with the changes predicted by the equation to provide a more complete 
prediction. 


2. Evaluation of Terms. 


The actual values used to evaluate the equation must be based on an 
adequate number of representative measurements and samples taken along the 
entire reach of shore for which the prediction is desired. Having carefully 
made these determinations, the engineer must then rely on the accuracy of the 
closure depth values given in the Appendix. These published values are based 
on the assumption that the profile closure depth is proportional to storm wave 
heights at the site. The proportionality constant has been estimated using 
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hindeasted, deepwater storm wave data of Resio and Vincent (1976) and actual 
profile changes measured in Lake Michigan (Hands, 1980). If the engineer 
feels that the published values seem too deep for a particular site (or per- 
haps too shallow) but lacks sufficient evidence to justify making a different 
choice, it would be easy and worthwhile to determine what effect the suspected 
error would have on the predicted response and on the conclusions that the 
prediction would support. 


An error in estimating closure will cause errors in both the numerator and 
denominator of the equation. In some cases, these two errors may reinforce 
each other; in other cases, they may compensate for one another, or they may 
even cancel and have no effect on the results or conclusionse Which situation 
exists can be determined readily if, as previously recommended, a number of 
shore-normal profiles have been obtained, and the following steps are takene 


First identify the point where the typical profile reaches the published 
closure depth, then the point where it reaches the alternately contemplated 
closure depth. Draw a line through these two points. If the line extends 
below the average height of the erodible backshore deposits, then overesti- 
mating closure depth will cause the equation to overpredict the response; 
underestimating closure will cause the equation to underpredict the response. 
If the line extends above the backshore deposits, overestimating closure will 
underpredict response and underestimating will overpredict response. These 
relationships will apply regardless of whether the response is a retreat or an 
advance of the shore. The possibilities are indicated in Figure 2. Inter- 
section of the line with the crest of the backshore deposits means that the 
choice between the two closure depths will have no effect on the calculated 
outcome. 


Predicted losses will be too small if the estimated 
UNDERCUTS DUNES closure depth is too small. 
le with 
WATER RISING Predicted losses will be too large if the estimated 
closure depth {s too large. 


Predicted gains will be too small if the estimated 
- closure depth is too small. 
with 
WATER EACCING Predicted gains will be too large if the estimated 
closure depth 1s too large. 


IF 
EXTENSION OF 
THE LINE 
THROUGH 
4) and d, 


OVERSHOOTS DUNES Predicted losses will be too large if the estimated 
fe. closure depth 4s too small. 


with 
WATER RISING 


Predicted losses will be too smal] if the estimated 
closure depth 1s too large. 


“Ach “ Predicted gains will be too large if the estimated 
losure depth {s too small. 
WATER FALLING © 


Predicted gains will be too small if the estimated 
closure depth 1s too large. 


Figure 2. Diagram for determining if a suspected error would weaken or 
strengthen arguments based on the sediment balance prediction. 
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VI. CONCLUSION 


A method has been presented for predicting one aspect to the long-term 
evolution of sandy shores. That aspect concerns displacement of the mean 
shore profile in response to long-term changes in water levels. This pre- 
diction method worked well in one hindcast instance involving a period of 
principally rising water (Hands, 1980). The degree of success found in future 
applications should be reported among Great Lakes engineers. 
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APPENDIX 


ESTIMATED DEPTH OF PROFILE CLOSURE 


Ay _ 
“es, 


The height of the fence shown above is proportional to the inferred 
closure depths obtained by relating extreme deepwater wave statistics and the 
results of repetitive profiling on the eastern shore of Lake Michigan (Hands, 
1980). Depths of profile closure beneath the stillwater surface are tabulated 


(in meters) on the following pages. 
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Profile closure depths (in meters) at Lake Superior study sites. 


SITE LOCATION LATITUNE LONGTITUNE CLOSURE 
HWEPTH [MI 

1 PIGEON FAY, tn “7.95 B9 U4 41.1 

2 GR4tat PORTAGE, MN 47.60 89.63 N@o® 

3 FARQUHAP KHNOK, tN 47,80 69,83 Ook 
4 BRULE RIVER, MN 47.680 90.08 12.4 

a GRAHL MAPATS, HN 4U7.66 90.28 12.4 
é POPLAR RIVER, MN 47.66 90.49 12.4 

7 CARLTON PEAK, MN 47.50 90.68 13.0 
8 TACONITE HARBOR, MN 47.50 90.91 12.6 

9 FAPTISM KIVER, MN 47. Su 91.451 12.6 
10 LITTLE TWO HARBORS, MN 47.20 91.32 12.6 
J1 AGATE KAY, ft 47,05 Pu oda 11.8 
12 ENIFE RIVER, MN 46.90 Palo Ck 11.6 
13 NULUTH, Mf 46.75 Pho PE VBenk 
14 RRULE POINT, WI 46.75 Palo Ck Poe 
13 TROW RIVER, WI 46.90 DA ed 8.2 
16 CRANKERRY RIWER, WI 46.91 91.30 8.8 
17 SISKIWIT KAY, WI 47.06 91.09 &.0 
18 POINT METOUR, WI 47.07 90.88 8.4 
19 ROCKY ISLANII, WI 47.08 90.67 8.6 
20 HARELE POTHT, WI 46.78 90.45 Fo tf 
21 SAXON HARBOR, WI 4U6.6u 90.45 8.2 
22 MORTREAL RIVER, MI 46.64 90.24 8.6 
23 PREGSHHE ISLE FIVER, HI 446.60 90. Ob 8.6 
24 PORCUPINE MOUNTAINS, fl 46,94 89.81 10.35 
20 PORCUPIME MONMTAINS, MI 46.9 89.41. 10.5 
26 CNTONAGON, iT 46.94 89.40 ial gal 
27 FOURTEEM Mie POINT, I 47.09 89.19 10.9 
28 atti) ISIE, tall 47.11 88.98 al oS 
29 REMRIGEE, tH) 7.24 83.77 11.3 
30 CALUMET, fil 47.36 82.55 alo a 
3] EAGLE RIVER, mI U7.5S 68. 34 10.3 
32 EAGLE HAREOR, I 47.53 88.14 DoF 
33 COPPER HARBOR, MI U7. 87.9) Dot 
34 SCHLATTER LAKE, MI U7.5S 87.70 °.0 
3S MONTTOU ISLANT, MI 47,36 87.49 11.3 
36 KEWEENAW POINT, MI 47.38 87.70 11.6 
37 POINT ISABELLE, MI 47, 2u EC. 94 TTS) 
38 TRAVERSE POINT, HI 47.0° 86.14 11.4 
39 PEQUAMING, HI 46.95 86, 35 PoP? 
40 HURON RIVER POINT, MI 46.95 87.91 Poe 
tI EIG kAY, AI 46.935 v.70 10.1 
u2 GARLIC POINT, MI 46,61 87.9 10.3 
43 MARQUETTE, MI 46.66 Coe? 10.7 
uy NEERTON, NI 46.66 87.08 10.9 

ral AU TRAIN FAY, MI 46.66 86.86 10.7 
US GRAND TSLAND, MI 46.66 66.635 stat ga 
uy GRANT! PORTAL POINT, MI 46.45 86.45 11.3 
ug AU SABLE POINT, MI 46.79 86,235 11.6 
4d GRANT! MARAIS, MI 46.78 84.02 11.6 
20 SUCKER RIVER, HI 46.78 85.80 11.64 
3 TEER PARK, MJ 46.78 85.60 11.8 
v2 LITTLE LAKE HARKOR, MI 46.77 B5.39 11.8 
oS CRISP POINT, MT 46.93 85.18 41.1 
oY PARADISE, HI 46.62 84.97 6.3 
Nike) POINT JROQUOTS, MI 4G6.U6 84.75 8.2 
34 GROS CAP, HI 46.61 Bu. 74 6.3 
a7 GOULAJS RIVER, MI 46.76 64.75 11.6 
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Profile closure depths (in meters) at Lake Michigan study sites. 


SITE LOCATION LATITUNE LONGITUDE CLOSURE 
MEPTH CMI 
1 STURGEON KAY, MI §5.64 85.32 ever 
2 HARKOR SPRING, HI 435.50 85.32 7.8 
3 FISHERMAN ISLAND, MI 45.37 &5.50 9.0 
4 TRAVERSE BAY, HI 45.23 85.72 8.6 
5 TRAVERSE CITY, HI 45.10 86.15 7.8 
é MANITOU, MI 4495 86.15 8.2 
7 PLATTE !.AKE, MI 44.80 86.17 9.7 
8 FRANKFORT, HI 44.68 86.39 10.9 
9 ARCADIA, HI 44.53 86.40 10.9 
10 ONEKAMA, HI 44,38 86.40 9.0 
ii MANISTEE, M1 GY, 24 B6.41 9.5 
12 RIG SAKLE POINT, MI &4y 09 86.63 9.0 
13 LUDINGTON, NI G3.94 86.64 sta 5 oh 
14 PENTWATER, MI 43.80 86.66 (0 4? 
15 LITTLE SAKLE POINT, HI 43.45 84.66 15 oS! 
16 RENONA, HI 43.52 86.47 ot oS 
17 MONTAGUE, HI Y3.36 86.48 12.0 
18 MUSKEGON, MI 43.2 P6.50 12.0 
iy GRANITE HAVEN, ‘HI 43.06 86.32 eae 
20 GRANDE RAPIUS, MI 42.93 86.33 11.8 
21 HOLLANI, MI 42.78 86.3 3U5t gat 
22 IlQUGLAS, HI Y2.64 86.35 SUSt 6 SS 
23 SOUTH HAVEN, MI 42.48 86.36 shat gph 
24 S SOUTH HAVEN, -MI 42,34 86.357 LO 
25 FKENTON HAREOR, MI 42.21 86.57 10.5 
26 S ST. JOSEPH, HI 42,06 86.58 10.5 
27 NEW KUFFALO, HMI 41.93 84.78 10.1 
28 MICHIGAN CITY, IND 41.78 86.99 10.9 
29 FBURNS HAREGR, INE 41.79 87.18 Lil gw 
30 CHICAGO, ILL 41.80 87.38 10.7 
31 CHICAGO SHIP CANAL, ILL 41.95 87.56 skal 5 XS 
2 EVANSTON, ILL 42.10 87.56 rato al 
33 HIGHLAND PARK, JILL 42.2 Bivens) LN oS 
34 WAUKEGAN, ILL 42.40 87.73 8.0 
35 KENOSHA, WI 42,54 87.73 10.9 
36 S RACINE, WI 42,69 87.71 8.2 
37 N RACINE, WI 42.83 87.70 sUOin, 
38 S MILUAUIKEE, WI 42.97 87.69 10.5 
39 MILWAUKEE, WI 43.12 &7.é68 10.1 
4&0 S PORT WASHINGTON, UWI 43.27 67.68 8.2 
41 PORT WASHINGTON, WI Y3.41 87.47 T of3} 
42 N PORT WASHINGTON, WI 43.55 87.646 8.0 
US S SHEHOYGAN, UI 43.69 87.45 8.4 
Gy N SHEROYGAN, WI 43.84 87.65 OFAS) 
4S HANITOWOC, UWI 43.98 87.64 8.0 
4S TWO. RIVERS, WI 44,13 87.43 Thal 
4? RAWLEY POINT, WI Qu, 27 87.42 9.5 
4B KEWAUNEE, WI 44,42 87.41 8.0 
49 ALGOMA, Ut GY .56 87.20 7.6 
50 STURGEON KAY CANAL, WI 4u,70 87.20 7.4 
Si JACKSONPORT, UI 44.84 87.18 9.9 
52 KAILEYS HARKOR, WI G4 .98 86.98 9.5 
33 N CANAL. LIGHT, WI 45.14 86.96 7.6 
S54 WASHINGTON ISLAND, WI 45.27 84.74 ors 
55 FISHERMAN SHOAL, WI 45.41 86.73 6.7 
56 ESCANALA, HI 4S.54 86.52 9.9 
3? POINT AUX KARQUES, HI 43.48 86.30 8.6 
58 N POINT AUX KARQUES, MI 45.83 86.29 6.7 
39 HANISTIQUE, HI 45.83 86.08 8.2 
60 PORT ISLAND, HI 45.81 65.88 7.8 
61 POINT PATTERSON, HI bS.94 R5.46 Si) 
62 MILLE COQUINS REEF, HMI usS.94 85.45 u.8 
463 SAULE STE MARIE, HI 45.93 85.25 5.0 
64 KREWOORT LAKE, HI 45.92 85.04 7.8 
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Profile closure depths (in meters) at Lake Huron study sites. 


SITE LOCATION LATITUDE LONGITUDE CLOSURE 
DEPTH £rI 
j PORT HURON, MI 43,02 Siena 13.0 
2 LAKEPORT, HI US.16 82.33 Ue & 
3 LEXTHWGTOGH, MI US. 31 82.335 Go 
4 PORT SANILAC, MI US .4uS” 62.52 UO oF 
ra FORESTVILLE, MI WA .59 He 9 DE 7.& 
6 HELENA, MI 43.74 82.52 Los 
i’, HAREGE EEACH, MI 43.88 Se 9 DL 12.0 
8 HURON CITY, NI 4.03 G2.71 Uo & 
9 PORT CRESCENT, MI uy, 05 63.11 12.4 
10 ENTRANCE SAGINAW KAY, MI uu 03 85.20 12.6 
iy TAWAS CITY, MI uu L17 63.3 13.0 
ee OSCOUS AU SABLE, MI 44,31 83.11 12.4 
AS GREENBUSH, I UU US SS) gdh dh 12.4 
14 HARRISVILLE: WI bu, 40 eS.1i1 13.0 
ja BLACK RIVER, mi tee, Puy fS.11 en tt 
14 QSSTHEKE, WI uu, 8S 83.11 12.8 
17 WORTH POINT, MI bo. OU 62.11 12.8 
18 ROCKPORT, I 45.18 63.50 12.6 
19 STONE PORT, MI uo. 42 63.20 WO oY 
20 ANAMS POINT, MI “S46 83.50 Po tt 
24 ROGERS CITY, MI US 4S £3.76 10.35 
22 HAMMOHT BAY, MI US. 61 65.90 2 oF 
23 CORDWOUT POINT, MI US. 7a 84.09 8.8 
24 POINT DOLOMITE, MI 4u5.89 84.09 11.3 
20 TETOUR REEF, MI bo. BS 85.90 10.5 
26 Wo ENT! DRUMMOND ISLAND, MI US. 75 63.70 10.3 
aul FALSE TETOUR CHANNEL, MT uo. 7o 83.50 11.3 
28 COCKBURN ISLANI, I 45.75 6835.30 11.8 
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Profile closure depths (in meters) at Lake Erie study sites. 


LOCATION 


MONKGE, MI 

CENAR POINT, OH 
LOCUST ROMmNi a Ok 
PORT Chunmoi Oli 
LAKESIME, OH 
HURON, Od 
VERMILIGHN, OH 
LORAIN, OH 

AVOH POINT, OH 
CLEVELANT, OH 

E GF CLEVELANT, @nH. 
FAIRPORT HAREOR, OH 
E GF FAIRPQRT HAREQR, OH 
GENEVA, QH 
ASHTARULA, OH 
CONNEAUT, OH 
GIRALM, PA 

BRE he 

f= OF (eRe, PA 

E GF NORTH EAST, PA 
WEST MEL 7 IN 

TAK IRE, NY 

ANGOLA, NY 
BUFFALO, NY 
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Profile closure depths (in meters) at Lake Ontario study sites. 


SITE LOCATION LATITUPE LONGITUNE CLOSURE 
DEPTH EMI 


J FORT NIAGARA, NY 4Y3.43 CP 50S &.0 
2 WILSON, WY US .uUS 78.83 8.6 
3 APPLETON, NY YS.uUS 78.63 Tod 
4 Winer (Il PPI 5 INN US LUS 723.46 Goal 
2 LAKESITIE PARK, NY US.US Wo oe 7.6 
6 KENTIIALL, NY 4S .47 Ca OF ine 
i HORTH HAALIN, NY US.4UT Ut BY 8.2 
8 WEST ROCHESTER, NY W3.48 Choe) 7.6 
9 EAST ROCHESTER, NY 43.37 C7 UG? 9.0 
10 PULTNEYVILLE, NY. WS. 47 ¢?.29 P oe 
1j SOMUS, NY 43.38 Chole Go 
We WOLCOTT, NY 43.358 76,87 wad 
13 FAIRHAVEH, HY YS,48 76.68 eo & 
14 OSWEGO, NY 43.63 76,50 11.8 
=f LACINA, RY u3.é3 76,28 Veok 
16 BELLEVILLE, RY uS.@7 0 Gis} 13.4 
17 GALLOO ISLAND, NY 43.91 76,990 10.1 
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